/iSR study of the quantum dynamics in the frustrated S = | kagome bilayers 
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We present /^SR experiments in the S = | kagome bilayer compound Ba2Sn2ZnGaio-7 P Cr7 P 022 
(BSZCGO(p)) and compare it to the isostructural SrCr9 P Gai2-9 P Oig (SCGO(p)), including for the 
latter new results for p > 0.89. Quantum-dynamical low energy magnetic excitations are evidenced 
in this novel compound. We study the evolution of the muon relaxation rate with p, T and field. A 
phenomenological model for the muon relaxation based on sporadic dynamics due to spin excitations 
in a singlet sea proposed by Uemura et al. [Phys. Rev. Lett. 73, 3306 (1994)]. is extended to all 
fields and T-range. Its connexion to the RVB picture is discussed, and we argue that such coherent 
states might mediate the interactions between "impurities" which induce the spin glass freezing. 

PACS numbers: 75.40.Gb, 75.50.Lk, 76.75.+i 



The stabilization of a spin liquid state, in which all 
kinds of magnetic correlation functions are short ranged, 
seems now granted theoretically in the S = | Heisenberg 
kagome lattice with nearest neighbour (NN) antiferro- 
magnetic interactions. The ground state (GS) is proba- 
bly of resonating valence bond type (RVB), built from a 
macroscopic number of singlet states [l| . If any, a gap be- 
tween the GS and the first triplet state is expected to be 
fairly small (~ J/20) in comparison with the exchange 
constant J. Most striking, due to both corner sharing ge- 
ometry and half integer spins, an exponential density of 
low lying singlet excitations, with energies smaller than 
the gap, ispredicted 0] and yields an extensive entropy 
at T = [4j. Finally, the nature of the first magnetic 
excitations, such as unconfined spinons, are still under 
debate 0. 

Among all corner sharing highly frustrated magnets, 
only a few are good candidates fulfilling the important 
conditions of the pure Heisenberg lattice with NN cou- 
plings. The combination of the weakness of the single- 
ion anisotropy and of a direct overlap exchange are cer- 
tainly the major advantages of the recently discovered 
chromate 5=1 kagome bilayer BSZCGO(p) @, and 
the long studied SCGO(p) pj, where p is the Cr 3+ net- 
work covering rate. Beyond the absence of ordering well 
below the Curie- Weiss temperature, the unusual large 
value of the specific heat unveils a high density of low ly- 
ing excitations 0, and its field independence suggests 
that the excited states are mostly singlets [jj. Moreover, 
their GS is found essentially fluctuating as proven by neu- 
tron experiments [l^, and /xSR for SCGO(p) [H [3, al- 
though an intrinsic spin glass (SG) signature is observed 
in susceptibility measurements pi [lj| . The origin of such 
a SG state in a disorder-free system still awaits for a 
complete understanding although recent theoretical ap- 
proaches catch some of the facets of this orig inal GS hj. 

In this context BSZCGO(p) appears as a unique com- 



pound : it has a more pronounced 2D character than 
SCGO(p), with a larger interbilayer distance H and all 
the magnetic sites lying on a corner sharing lattice [22j |. 
In addition, a twice lower SG temperature T g has been 
reported Q and a new type of defects other than non- 
magnetic substitutions (likely bond defects), dominates 
the low-T susceptibility jj . This offers the possibility of 
probing the influence of T g and the relevance of dilution 
on the evolution of the T — > dynamics, by a simple 
comparison of two quasi-identical kagome bilayers. 

In this Letter, we present the first Muon Spin Relax- 
ation (/xSR) study of the spin dynamics in BSZCGO(p) 
and revisit the case of SCGO(p) for low dilutions now 
available. This technique has proven to be a front tool 
for the study of quantum dynamical states in a vast fam- 
ily of singlet GS systems. Whereas specific heat is sen- 
sitive to all kinds of excitations, including likely domi- 
nant low lying singlets at low-T muons (/i + , S = g) 
probe magnetic excitations only. Their relaxation rate 
A is linked to the spin-spin time correlation function 
A = 1/Ti - / o °°(S(0)-S(i)) cos{ ltl H LF t)dt, where H LF is 
the external field applied along the /i + initial polarization 
and 7 M is the /x + gyromagnetic ratio. The typical time 
window (10 ns-10 [is) allows to observe spin fluctuations 
which set in the specific range of frequencies (~ 10 9 Hz) 
in between NMR and neutron experiments. Not only in 
SCGO [I1E3 and in other kagome systems [HE3 but 
also in other spin-singlet compounds [l7j . the signature 
of the quantum dynamics when T — > could be identified 
to the occurence of a plateau of A and a field-independent 
initial Gaussian shape of the /i + polarisation P z (t). Ex- 
tending the original guess of Uemura et al. to explain 
this exotic fluctuating GS we propose for the first 
time a phenomenological model which explains both T- 
and field-dependence of P z (t) for all p in these kagome 
bilayers. 

In order to illustrate the evolution of the properties 
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FIG. 1: a,b) T-dependence of the polarization P z (t) in 
BSZCGO under weak (10 G for p = 0.43 and 100 G for p = 
0.97) longitudinal field. Notice the different time scales for 
both samples. c,d) ffiF-dependence at base-T. The line for 
p = 0.43, ZF, is a square-root exponential multiplied by the 
Kubo-Toyabe function in zero field (Ajvd ~ 0.08/is" 1 ), to 
take into acount the nuclear dipole contribution. The other 
lines are fits described in the text. 



with dilution, we first focus on the qualitative differences 
(Fig. ^) between two typical low and high occupations p 
of the kagome bilayers. Starting from high-T (T > 6 K 
in Fig.^a,b)), a conventional paramagnetic behavior is 
found for both cases, with an exponential variation of 
P z (t) in the dense case and square root exponential in the 
dilute case. The evolution of P z (t) at low-T is markedly 
different. For p = 0.97, the relaxation rate increases 
by more than two orders of magnitude to reach a T- 
independent value for T < T g w 1.5 K 5| (Fig. Eft))). 
For T = 0.03 K, P z (t) displays a shape in between ex- 
ponential and Gaussian at early times and reaches zero 
value at long times (Fig. ^d)). It also displays a much 
weaker field dependence than for a frozen magnetic state. 
Indeed, for A ~ 2 ^s" 1 , a 500 G field would completely 
decouple the n + polarisation E3 . Following: the argu- 
ments initially developed for SCGO(p) all these 
properties are the first indications of a quantum dynam- 
ical magnetic state in BSZCGO(0.97). For p = 0.43, we 
also find a dynamical state but, at variance with the pre- 
vious case, only a weak T-dependence is observed. Also, 
the polarization displays a square root exponential decay 
for all T (Fig. da)) and for any field H LF (Fig. [He)). A 
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FIG. 2: Comparison of the evolution of the n + relaxation 
rate A vs. T in BSZCGO(0.97) and SCGO(0.95), from a 1/e 
analysis. Inset : recovery of a part of the asymmetry at low-T 
and long times in SCGO below T g (log- log scale). 



weak plateau of the relaxation rate is still observed be- 
low 0.5 K. However the relaxation rate increases weakly 
but continuously when T decreases for p = 0.3, with a 
very weak Hpp -dependence. This is typical of the pure 
paramagnetic fast fluctuations limit for diluted magnetic 
systems, although strongly frustrated antiferromagnetic 
interactions are still reported [f|. 

In a first step of analysis, we estimate the fi + relax- 
ation rate A using the \ point of P z (t), as discussed in 
[l2|. This allows to single out information about two 
important issues using a simple model-independent anal- 
ysis : the role of the SG-like transition and the impact 
of the non-magnetic dilutions on the dynamics. 

From high-T, A increases by more than two orders 
of magnitude down to T g for both purest samples, to 
reach the relaxation rate plateau Xt^o for T < T g . 
With a T-scale twice larger for SCGO (0.95) than in 
BSZCGO(0.97)^\(T) perfectly scales on the T axis as 
shown in Fig. [21 23] . This ratio is very close to their T g ra- 
tio of 2.3(2), which points to a link between the formation 
of the SG-like state and the presence of quantum fluctua- 
tions. This seems to be a quite common feature of various 
systems with singlet GS systems [HElElElEzl 

In Fig. [3| we report the variation of At— >o (T) for vari- 
ous BSZCGO (p) samples and compare them to SCGO(p) 
[ill . Although additionnal p-independent defects are 
dominant in BSZCGO (p) and evidenced from a large 
Curie- like term in macroscopic susceptibility at low-T [fj| , 
it is very surprising to find a very similar quantitative 
At^o variation in the whole p range for both systems. 
Therefore we conclude that only defects related to dilu- 
tion influence the relaxation rate. 

In a classical framework, a coplanar arrangement 
with zero-energy excitation modes, involving spins on 
hexagons, was proposed in the literature to be selected at 
low-T 18]. Inelastic neutron scattering experiments on 
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FIG. 3: p-dependence of \t^,o f° r BSZCGO (open symbols) 
and SCGO (full symbols). Continuous (dashed) lines are ~ 
p 3 + 1.9(3)p 6 (cxp 3 , HI) fits. 



the spinel compound ZnCi'204 are the first recent indica- 
tion in favor of such excitations |19|. Therefore, one ex- 
pects the fi + relaxation to scale with the number of fully 
occupied hexagons, cx p 6 . The BSZCGO(p) data alto- 
gether with our SCGO0.89 < p < 0.95 samples indicate 
that Xt^o (p) is well accounted for by adding a dominant 
p 6 term to the cx p 3 term proposed for SCGO(p < 0.89) 
in 0] (Fig. |3J ■ This certainly indicates that the relax- 
ation is not induced by single spin excitations but rather 
collective excitation processes extending at least on tri- 
angles and hexagons. 

One of the pending challenges in the analysis of /iSR in 
kagome frustrated antiferromagnets is to find a theoret- 
ical model reproducing P z {t) for all fields and all tem- 
peratures. In the quantum framework of a singlet GS, 
one needs unpaired spins in order to generate magnetic 
excitations responsible for the /i + relaxation. Such ex- 
citations can be ascribed by spin | unconfined spinons. 
Their location vary in time with no loss of coherence of 
the excited state. Hence, a given /i+ couples to a spin 
only a short fraction ft of the time t after implantation 
and one can use a model of "sporadic" field fluctuations 
to describe P z (t) . This was the initial guess by Uemura et 
al. to account for the weakness of the field dependence of 
P z (t) at low-T in SCGO(0.89) [H|. This polarization is 
given by a "sporadic" dynamical Kubo-Toyabe function 
which can be approximated by the Keren analytical func- 
tion P?(ft, A, lhy. v) = Pf{t, /A, fH LF , fv), valid in 
the limit v > A 20], quite appropriate here. A/j^ is the 
local field created on the /i + site by a paramagnetic neigh- 
bouring released spin, with a corresponding exponential 
time correlation function, exp(— ist), characterized by a 
fluctuation frequency v. A is related to the /i + location 
in the unit cell and its value is therefore p-independent 
for a given system but is expected to vary from BSZCGO 
to SCGO. The Gaussian at early times, the weakness of 
the field dependence and a dynamical relaxation down to 
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FIG. 4: Fitting parameters of Eq. a,b) Hlf- and 
T-dependence of x for BSZCGO(p) (open symbols) and 
SCGO(p) (full symbols). c) p-dependence of /. d) In- 
dependence of A'. 



are altogether related to the / factor and v ~ A. In ad- 
dition to this sporadic relaxation, we found that a more 
conventional Markovian relaxation channel needs to be 
introduced in order to fit the long times tail (t > 3 /is) 
for all fields and T, which yields an exponential variation 
of P z (i). We therefore write 



P t {t) = *P*(t, /A, fH LF , fv) + (1 - x)e 



-A't 



(1) 



where x represents the weight of the short time sporadic 
dynamical function, associated with the spinons dynam- 
ics. In the following, we first detail how all these param- 
eters can be reliably deduced from the data and sketch a 
physical picture consistent with our results. 

We first present our analysis for T <C T g in BSZ- 
CGO^). In order to limit the number of free parameters, 
we make the minimal assumption that the external field 
does not influence the dynamics of the coherent spinon 
term, v and A values are shared for all p and Hlf and 
are set by our data taken at low fields {Hlf < 500 G) 
for various p. f is adjusted for each p and its variation 
accounts for the evolution of Xt^q{p). For low fields, x 
is found close to 1 for the purest samples (Fig. IHa,b)), 
making A' a non-relevant fitting parameter. On the con- 
trary, the high field data enable us to track the evolution 
of x with Hlf and to determine a value of A'. We could 
not extend the fits below p = 0.71 since the weak field 
dependence prevents an unambiguous determination of 
the parameters. 

We find perfect fits of our data (Fig. IHd)) with v ~ 
1000 /is -1 , A ~ 350 /is -1 ~ 7^ X 4 kG and an average 
value of / ~ 0.006. We find a nearly flat p-independent 
long time relaxation rate for T <C T g (A' ~ 0.05 ^s _1 , 
-A- in Fig. E^d)). An important finding is that x is of 
the order of p at low fields (Fig.^Ja,b)). This may be re- 
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lated to theoretical computations showing that the corre- 
lations are enhanced around spin vacancies in the kagome 
lattice [U, which would destroy locally the liquid state. 
Besides, x decreases appreciably for H^p ~ 10 kG what- 
ever the value of p. We can associate this decrease to 
the existence of an energy scale ~ 1 K, which is of the 
order of T g . Finally, we observe a linear variation of / 
with p (Fig. Efc)) and / tends to vanish around p ~ 0.5, 
a limit consistent with our classical approach (Fig. |2J. 
This indicates that even far from the substituted sites, 
the coherent state is somehow affected, e.g the density of 
spinons could be smaller. 

As suggested by the similar p-variation of Xt^o in 
both BSZCGO and SCGO, we assume that the excita- 
tion modes are identical in both systems, i.e., / and v are 
kept the same for comparable p. We find A ~ 1200 /xs _1 
for our SCGO samples (p > 0.89), in agreement with pre- 
vious work on SCGO(0.89) |ll|. In order to cross-check 
these values of A, we use the high-T (~ 50 K) data for 
BSZCGO (SCGO), A - 0.01 ^s" 1 (0.03 ^s" 1 ). In the ap- 
propriate paramagnetic limit, the fluctuation rate can be 
estimated from v ~ yfzJS/kBfo ~2x 10 13 s" 1 11] using 
a coupling J = 40 K [(| and an average number z = 5.14 
of Cr NN 0. We extract A = y?Jvj2 ~ 300 /is" 1 
(600 ^s^ 1 ), consistent with our fits. It is quite rewarding 
to find that a common physical picture underlies all the 
sets of data at low-T for both kagome bilayers. 

We now extend the fits of P z (t) to the whole T-range 
(Fig.^b)), fixing / to its low-T value in order to limit the 
number of free parameters. As expected from the change 
of shape around T g [ill flij ] , the weight x of the sporadic 
term decreases to finally enter a high-T regime (Fig.^b)) 
where the p, + relaxation becomes exponential (x — » 0). 
Notice the similarity between the field and T-dependence 
of x which indicates that the sporadic regime is destroyed 
when the energy is of the order of T g . Fig.Efb) displays 
a sharp crossover from one state to the other, between T g 
and ~ 3T g , corresponding to the A steep decrease (Fig.|5J, 
for both BSZCGO(p) and SCGO(p). For T > T g , v and 
A' (Fig. Ud)) decrease by one order of magnitude up to 

At high-T, it is natural to think in terms of paramag- 
netic fluctuating spins. Lowering T, A' seems to diverge 
at T g and below, the weight of the exponential term at 
low fields (1 — x ~ 1 — p) could correspond to localized 
frozen spins, reflecting a glassy component. It is note- 
worthy that T g does not increase with 1 — p but rather 
decreases, at the opposite of the case of canonical spin 
glasses. We could further confirm the existence of such 
a frozen component in SCGO(p = 0.95 — 0.89), since a 
clear recovery of a small part (~ 4%) of the asymmetry 
is found at long times (e.g. Fig. EJinset)). It is observed 
for t ~ 100 A -1 , which is out of the experimental range 
in the case of BSZCGO (p). 

To summarize, the picture based on a coherent RVB 
state, which magnetic excitations are mobile fluctuating 



spins 5 on the kagome lattice, explains well the data, pro- 
vided that (i) these excitations can be generated even for 
T — > 0. This underlines the smallness of the "magnetic" 
gap, if any, typically A < J/1000 ; (ii) an energy scale re- 
lated to T g (fields of the order of 10 kG or, equivalently, 
temperatures of the order of T g , which vary very little 
with p) are high enough to destroy the coherent RVB 
type state ; (iii) the substitution defects are accounted 
for by an additional classical relaxation process. 

The absence of energy scale for spin I excitations could 
explain quite well why the transition to the SG state 
is fairly independent of p. Indeed, in this framework, 
spinons could mediate the interactions between magnetic 
defects localized around spin vacancies. This is corrobo- 
rated by all the scaling properties in T/T g . The transi- 
tion to the SG state corresponds then to the formation 
of the coherent singlet state rather than any interaction 
strength between defects. 
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